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Abstract 

The quaternion equation X n = A is solved for any integer number 
n > 2. A is a given quaternion with komplex numbers as its elements. 
We use the isomorphism between quaternions and (4, 4)-matrices to solve 
this equation. 



1 Notations 

latin lower case letter - complex number 
i denotes here the imaginary unit (i 2 = —1) 
bold latin lower case letter - a 3-dim. vector, 

0-2 

a 3 

bold latin capital letter - a quaternion, 
( ao \ 

a 2 \ a 

V 0? / 

A quaternion is often written as clq 
the above definition. 

-ai 
ax a Q 
a 2 a 3 
\ a 3 -a 2 

The inner product - ab := a\b\ + a 2 b 2 + 0363. 



i. e. A = 



a x i +a 2 j 
( ao 



a 3 k but in this article I use 



capital Gothic type - a matrix i.e. 21 = 



-a 2 
-a 3 
a 



-0,3 \ 
a 2 

ao ) 



The quadrat of a vector - a 2 := a\ 
The quaternion product - 



AB = 



/ a ° \ 




fb \ 




I a b Q 


- aih 


- a 2 b 2 


- a 3 b 3 \ 






bi 




a h 


+ aib 


+ a 2 b 3 


- a 3 b 2 


a 2 




b 2 




a b 2 


- aib 3 


+ a 2 bo 


+ a 3 bi 


V 03 J 




\ b * J 




\ aob 3 


+ aib 2 


- a 2 bi 


+ a 3 b Q J 



or 



AB 



«o 
a 



b 
b 



aobo — ab 
aoh + 60a + axb 



(1) 



(2) 



1 



The isomorphism between quaternions and 4x4 matrices // : Q — > M 4x4 is: 



( a ° \ 




( ao 


-ai 


—a-2 


-«3 \ 


ai 


) = 


a\ 


ao 


-a 3 


0-2 


0,2 






a 3 


ao 




V «3 J 




\ a 3 


-a 2 


a i 


«o / 



There are another isomorphisms too but they do not yield another result. 

The characteristic polynomial of a matrix: c(y) :— det(y<£ — 21). 

The minimal polynomial of a matrix: m(y) :— annihilating polynomial of least 

degree. 

The adjugate to (y<£ - 21): := adj(y<£ - 21). 

This relation is valid - (y£ — 2l)2J(y) = c(y)€ 

The reduced adjugate to (y€ — 21): €(y) with (y€ — 2l)£(y) = m(y)<£ 



2 Method of solution 

We have to solve the quaternion equation 

X" = J 



(3) 



with n > 2, n e N. The way to the solution consists of two steps: 

1. Transform the quaternion equation into the isomorphic matrix equation. 



Instead of the equation X™ = A mit A = 



equation 



with 21 = 



f a ° \ 

ai 

02 
V 03 J 



ao 
a 



X" =21 



we solve the 



(4) 



V 



ao 


-ai 


a i 


a 


«2 


«3 




-02 



-a 2 
-a 3 
a 
ai 



-«3 
a 2 

-ai 
«o 



The solution set is all matrices X, satisfying this equation. 
2. From the solution set {X} we select the subset 
Xq = {X\\ X isomorphic to a quaternion}. 

The isomorphic relation provides the solution set of the quaternion equation 
directly. 

If we want to solve equation |(4]) , then we only have to transcribe the solution 
from a good textbook about the theory of matrices. Therefore we regard at 
first the structur of the matrices which are isomorphic to quaternion (we call 
they q-matrices in the following). This preliminary work is necessary, because 
we need the Jordan normal form of the q-matrix 21. 



3 Jordan normal form of q-matrices 

All of the occuring numbers ai, Xi etc. are komplex numbers (elements of C). 
Solving problem: Compute the transformation matrix il and the Jordan matrix 
Z in the equation 21 — itjil -1 for a given q-matrix 21. Q-matrices are of the 
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form: 



21 : 



ai 
\ a 3 



a 

«3 

-a 2 



-a-2 
-a 3 
a 
a i 



-a-3 \ 
a 2 

ao J 



The characteristic polynomial of the q-matrix 21 is c(y) = (y 2 — 2aoy- 



The minimal polynomial of 21 is m(y) = y 2 
adjugate matrix 

/ V - a -ai -a 2 -a 3 \ 
a>\ V - a -a 3 a 2 
a 2 a 3 y - a -a\ 
\ a 3 -a 2 ai y - a J 

The solutions of the characteristic polynomial are 

j2 ~. — llt — n „ — <j\/ q2 



-a 2 ) 5 



2aoy + Oq + a 2 and the reduced 



(5) 



(6) 



21 = 3 = 



«o 



/ 1 \ 
10 
00 10 
y 000 1 J 



(7) 



yi = y 2 = ao + iVa", y 3 = y 4 = a 

There are three cases for the Jordan normal form of 21 = itjll -1 : 

1. a = 0. 

In this case 21 is of the representation 

/ a Q \ 
a 
a 
\ a / 

The matrix 21 already has Jordan normal form with (after reducing) the minimal 
polynomial m{y) = y — ao and the reduced adjugate matrix 

/ 1 000 \ 
10 
10 

\ o o o i / 

2. a ^ und a 2 ^ 

Here the minimal polynomial has 2 simple roots. The corresponding Jordan 
matrix is diagonalizable: 



£(y) = € = 



I ao 



■6 = 



\ 



a 






a - n 








a - i\ 



\ 



J 



(«) 



The transformation matrix il is assembled by the columns of the reduced ad- 
jugate matrix £(y) (JH). Because a ^ and a 2 ^ 0, we have to differentiate 
between three cases: 
a) a 2 + a\ ^ 



il is assembled by the first both columns of €(ao + iva 2 ) and the last both 



columns of €(ao — iva 2 ) Then we have: 



«3 



-<Zi 

ix/a 2 
a-.i 

-a 2 



-a 2 
-0,3 
—iva 2 
ai 



-a 3 

«2 
-<Zi 

-iva 2 



(9) 
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and 



it- 1 = 



1 

2(a 2 +0 2 ) 



/ -i(a\ + 4) 


2 



V -a 3 VsP 



a2Va- 
— i(a2 + a|) a^Va 2 

ia\a,2 02 Va 2 + 20103 



2 — iaid3 a3\/a^ + «aia2 \ 
iaia2 — a2V / a 2 " + 20103 



a3\/a 2 +«aia2 *(« 2 



2 

2 t- "3 




4) 



o 

*(«! + a i) 



b) a? 



it results from the first and third column of £(a + iva 2 ) and the second and 
fourth column of £(oq — zVa 2 ) Then we have: 



It: 



( iva 2 
ai 
a 2 
o 3 



-o 2 - fl 3 

2 



-o 3 
iVa 2 

Oi 



and 



it- 1 = 



-ai 



1 



«2 

-ai 
-iVa 2 



(10) 



-ai 



i(a 2 + 4) 


,/^2 



2(a? + o§), 
ai%/a 2 + ia 2 a 3 



?aia2 



ia 2 a 3 +i{a\ + 4) 
ia\Ch2 



-z(aj + 4) 

>2 



a 3 va^ + ^ala2 
-aiVa 2 + 2a 2 a 3 



asVa 2 - zaia 2 \ 
aw/a? + ia 2 a 3 


i{al + 4) 



J 



V -asA/a 2 - 
c) a\ + a\ ^ 

it results from the first and fourth column of €(a + zVa 2 ) and the second and 
third column of £(a — iVa 2 ) Then we have: 



it: 



f Wa 2 -a 3 

dl 0,2 

a 2 —a\ 



\ 



a 3 



and 



it" 1 = 



— di —02 

-iVa 2 —03 
03 —iva 2 
— a 2 ai 

1 



(11) 



V 



-or 
-o 2 - 



-i(a? + a|) 


/ a 2 — ia 2 a 3 z(of + a 2 ) 



2(a 2 + a 2 )Va 2 

12 



aiVa- 4 — ia,2a 3 0,2V a* + iaia 3 
a 2 Va 2 + iaia3 — aiVa 2 + 40203 





-i(fli + a|) 



\ 



-a 2 va z 



«Oi03 



'a 2 + zoi03 







aiV / a 2 +«a 2 03 / 





i{a\ + 4) 
3. a ^ und a 2 = 

Here the minimal polynomial is m(y) = y 2 — la^y + 4 — {y ~ a o) 2 - The only 
(of multiplicity 2) solution is y\ = ao- The corresponding Jordan matrix of 21 is 
then: 

( a 1 \ 
a 
a 1 
\ s / 



3 = 



(12) 



Autor: Jochen Hans 



4 



We use the reduced adjugate matrix €(y) and the derivative: 

-as \ 



( y-a -ai -a 2 

ai V - a -o 3 a 2 
a 2 a 3 y - a -a\ 
a 3 -a 2 ai y - a J 



V 

/ 1 \ 
10 
00 10 
\ 000 1 J 

to compute the transformation matrix il. We have to differentiate between three 

cases too: 

a) a\ + al ^ 

il results from the first and the second column of £(oo) and of € '(oo). 



il 



/ 1 -oi \ 

ai 1 

a 2 a 3 

\ a 3 -o 2 J 



(13) 



and 



iT 



V 



a 2 a 3 \ 

aia 3 —a\a 2 

a 3 — «2 

a| + a 3 — aia 2 — ftia 3 y 



b) af + a\ ^ 

il results from the first and the third column of £(ao) and of € '(do). 



il = 



and 



il~ 



^0 



1 



a{ + a 3 



/ 1 

ai 

a 2 
\ a 3 



a\ + a\ 





-02 \ 

-as 

1 

01 o y 

Oi 

-a 2 a 3 
-a 3 

-aia 2 a\ 



(14) 



03 

aia 2 

Oi 

-a 2 a 3 



c) a? ■ 

il results from the first and the fourth column of <£(ao) and of € '(do). 



il: 



/ 1 -o 3 \ 

ai a 2 
a 2 — ai 

V a 3 o o i y 



(15) 



and 







ai 



il- 1 = 



a\ + a\ 



a\ + a 2 a 2 a 3 



a 2 
-aia 3 



a 2 — ai 
y -aia 3 -a 2 a 3 a\ + a 2 J 
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Why do I compute the transformation matrix il and the Jordan matrix 3 so 
extensive? In order to show that it and Z are uniquely determined by the 
four components of the quaternion A. Futhermore we found that a q-matrix is 
similar to one of three Jordan matrices Q , ® or (fT2| . Now we can go about 
the solution of (0]). For it we remove the appropriate chapter from a suitable 
textbook. 



4 Solution 

We also consult pQ chapter 8.6 at page 243. There we read that we have to 
make yet another case differentiation, the differentiation between invertible and 
singular 21. Now let us begin with 



4.1 Invertible matrices/quaternions 

1. Case - The matrix 21 is of type $7§. 
It can be written as 



21 = a 



(16) 



/ 1 00 \ 
10 
10 
\ 00 1 / 

Furthermore we need all invertible matrices 53 (chapter 8.6 from p]) which are 
commutable with 21 also satisfy the equation 2193 = 5321. Because 21 = ag(B the 
above equation holds for any invertible matrix 23. Let Wi be an arbitrary n-th 
root of qq. The general solution of ((!]) is then: 



X = 93 



\ 



/ w x 
w 2 
w 3 
\ Wi J 



ar 



(17) 



Here only the solution matrices are of interest, which are isomorphic to a quater- 
nion. Than there are the the following cases: 
a) wi = W2 = W3 — W4: Then is X = 23wi(£93 _1 = wi*£. 

(I) 

Solution quaternions X are any X = wi ^ . There are n solutions accord- 

v / 

/oq (This is the trivial solution.). 



ing to the n different values of 
b) w\ = W2 und W3 = W4 

Here w\ and w 3 are two different n-th roots of ag, t. i. w™ = 103 = ao but 
w\ 7^ W3. So that X is isomorphic to a quaternion, we need to construct il 
according to |(9|) , (fT0|) or ifTTj) . If we select any numbers for xq — X4 then we can 
compute the il without restrictions. The general solution is 



X = 23 



V 



wi 

W! 

w 3 




\ 



93 



— 1 



(18) 



w 3 j 
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with arbitrary invertible matrix QJ. Therewith X is isomorphic to a quaternion 
(of type JSl) we have the conditions 



Xq + iv x 



Wi 



U>3 = Xq — IV X 



(19) 



and the matrix QJ has to be an U-matrix according to , (fit)]) or (fTTjl . Because 
w\ and W3 are ^ we always can compute the corresponding u-matrix and 
X0...X4 only need to satisfy lfT9|) . The corresponding quaternion set {X} results 
into 

{X} = {X I x = hw x + 103) A x 2 = x\ + x\ + x\ = ~(W! - w 3 ) 2 } (20) 

There are (Zj of this solution sets (accordingly to the numbers of the possible 
combinations of w\ and w 3 ). 



1. Example: n=2, to solve the equation X 2 = oqE = 

1 



ftfi 




The simple 



solution is X = i-v/ao 







The second solution (according to b))is of more 



interest. The two different roots w\ and w 3 are: wi = ^/oq and W3 = —^/oq. So 
that the solution set is {X} with xq = and x 2 



-<Zq 



2. Example: n=3, to solve the equation X J = aoE = 



a 




The trivial 



solution is obviously X = tfao E, whereas for ^/ao tree values are possible. The 
solutions according to b) are ( 2 ) = 3 different sets, which are computed in the 
following. The 3 different 3-rd roots from ao are: Wi = tfao, u>2 = w\ (— 1+| v3) 



and W3 = wi(- 



;VS)- Then the combinations result: (wi,w 2 ), (wi,w 3 ) and 



(w2,w 3 ). This provides the following three solution sets: 



H(-f 



|V3) 



(wi,w 2 ) ^> so = wi(| + |V3) x 2 
(wi, w 3 ) =>• .t = - |V3) x 2 
(w 2 ,w 3 ) =>• x = ~\wi x 2 = 

2. Case: The matrix 21 is of type ([8]). 

Using a transformation matrix of kind ([9]), (fT0|) or l[TTj) the 21 can be written as: 
21 =H3ir 1 . 

2 is a diagonal matrix. An arbitrary, with 3 commutable, invertible matrix 53 
is of kind: (q.v. [I], chapter 8.2) 



( vi v 2 ^ 

W3 W4 

v 5 v e 
\ 0d 7 « 8 / 



(21) 



Then the general solution of (@]) is (q.v. [TJ, S. 245) 
/ voo + i\/a 2 



X = U2J 



ao + iv a z 






V 7 



ao + iv & 





ao — iva z 




ST 1 !!" 1 



a 



(22) 
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In this general formula different values can be inserted for the n-th root in the 
1. and 2. as soon as in the 3. and 4. row (accordingly the multiplicity of 
the n-th root). We need only X-matrices which are isomorphic to a quaternion. 
Therefore roots of the 1. and 2. row as soon as the 3. and 4. row of 3 have to 
be equal respectively. But then 



■6x — 



V 



a + i\ 










\/ao +W a 2 
Vflo-i 




\ 



(23) 



is commutable with 2J and the equation 3a 
/ Voo + iVa 2 



X = il 



V 





y/a + i\ 





Voo - iVe? J 
is true. X is 

\ 



\/a - n 




(24) 



\/a - iVa 2 " / 



The transformation matrix il is as same as that of 21. Then we can write X as 
quaternion as follows: 



2 

J__ 

"2* v 



X=\t+iV? +^iV^| ,- fl ). (25) 

n values can be inserted for the n-th roots so that follows altogether n 2 solutions, 
(every value of \/ao + W a 2 can be combined with every value of \/ao — iV a 2 .) 
here an excample: 



Let A 



76 

24i 
12i 

v -4i* y 



n = 3, to solve the equation X 3 = A. The two eigenvalues 



of 21 are with Va 2 = A9i: 

a a + iVt£ = 27 and a — iVs^ 



3 

2' vu ' 2 2 

written down as: 



§2\/3 as soon as 5; 



125. The cubic roots have the values: 3; — | + 
| — |i\/3- The 9 solutions are 



■f + fn/3; 



Xi 



7 


-12i 




/ 


\ \ 

Yli 






4 \ 

24i 




49 


+ 5 




49 






49 




-6i 




6i 






Yli 




49 






49 






49 


V 


41i 

98 / 






-41i J 
98 / 






-41i , 
49 / 



(-1+ 



fiV5) 



i 

2 

-12i 

49 
-6i 

49 
41i 



+ 5 



_L 

2 
12i 
49 
6i 
49 

-41j 



fiV5) 



i 

2 

-12i 
49 



49 
41i 



_L 

2 

12i 
49 
6i 
49 
-41j 



7+3z\/3" 
4 

1873+78t 
49 

9\/3+39a 
49 

-123v^3-533t 
196 

7-3^73 
4 

-1873+78t 
49 

-9^+398 
49 

123\Z3-533t 
196 
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1 




2 


\ 


-12i 




49 




-6i 




49 




41t 




98 



\ 98 



+ (- 



- 5 _i_ « , 
2 T 2* 



V3) 





1 ^ 

12i 






49 






6i 






49 




V 


-41? ( 
98 / 





( * \ 



x 5 = (- 



l 

2 

-12/ 

49 
-6/ 

49 
41/ 



2 

-12/ 

49 
-6/ 

49 
41/ 



+ (" 



X 7 = 3 



1 




2 


\ 


-12/ 




49 




-6i 




49 




41j 




98 



i-s/3) 



1+5/V3 
4 

-30\/3-78i 
49 

-15\/3-33/ 
49 

205\/3+451/ 
196 



_L 

2 
12i 
49 
6i 
49 
-41/ 
98 

_L 

2 
12/ 
49 

ill 

49 
-41/ 
98 



2iv / 3 \ 



-12\/3-12i 
49 



49 

41V3+41/ 
98 



*V3 \ 



-48\/3-12/ 
49 

-24V3-6/ 
49 

164y / 3+41/ 
98 



l~5iv/3 



\ Mi 

\ 98 



+ (-§ - |*V3) 



2 

12i 
49 

ill 

49 
-41/ 
98 



\ 98 / 



49 

15\/3— 33/ 
49 

-205%/3+451/ 
196 



Xo 



i 

2 

-12/ 

49 
-6/ 

49 
41/ 



1 

2 

-12/ 
49 
-6/ 
49 
41/ 



»V3) 



+ (- 



»V3) 



2 
12/ 
49 

111 

49 
-41/ 
98 

_L 

2 
12) 
49 
6i 
49 
-41/ / 
98 / 



\ 



V 



48V3-12/ 
49 

24V3-6/ 
49 

-164\/3+41i 

98 / 

-2-2n/3 \ 

12V3-12/ 
49 
6^-6/ 
49 

-41\/3+41/ 
98 



/ 



3. Case: The matrix 21 is of type l|T2j) . 

Here the transformation matrix is of kind (fT3|) , or lfT5l) and 3 is of kind ffl2|) . 
We have the equation 21 = il^il -1 . An arbitrary, invertible, with 3 commutable 
matrix has the form (according to [T], Kapitel 8.2): 

^ ui v 2 v 3 V4 \ 

Vi V3 

W 5 W 6 W 7 W 8 

\ t> 5 v 7 ) 
The general solution of (J4j) is then (with the abbreviation w 

( - 1 



27 = 



(26) 



ao) 



f'i 







wi 
1 



\ 



w 2 



V 



(27) 



/ 



w 2 

We need only the solutions X, which are isomorphic to a quaternion. In this 
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case the two roots u>i and W2 have to be equal. It follows the solution matrix: 

(28) 



X = SW 



( W —krrj \ 

w 

W —^T 

\ w 



Or 1 !! -1 



In this formula the w-matrix is commutable with 03 too and X reduces to 

/ to —kr^r \ 
to 
to — Itt 

\ to / 

The retranslation into the quaternion notation provides then 



X = il 



it" 1 



(29) 



X = 




10 
na 



( \ 

V «3 J 



(30) 



Thereby we have to = ^/Eq and there are n different values possible. 
( 1 \ 

Example: Let A = ^ , n = 4. Then it follows to = \/T with the solution 
V 5 / 

set {1; —1; i; —i}. The 4 solutions of X 4 — A are: 





/ 






( 


-1 
-H 


\ 


>-l — 




i 


x 2 = 




— i 






V 


I J 




\ 


5 
4 


) 




( 






( 


— i 








3 






3 


3 = 




4 
-1 


x 4 = 




4 
1 






\ 


¥ j 




\ 




) 



4.2 Singular matrices/quaternions 

Singular quaternions can too be divided into the three classes 0, JH]) and (fT2|l . 
For 21 we can make the same case differentiation as above at the invertible 
quaternions. Then we get different solutions as in the case of invertible quater- 
nions. 

1. a = 0. (Type ©) 

So that A, resp. 21 is the Zero quaternion, resp. the Zero matrix. The following 
equation is to solve: 

/ ^ 

„ = 0000 


\ J 
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Evidently the Zero matrix is commutable with any matrix. The solution of this 
equation is of the general form X — T-O^T -1 with a transformation matrix 1, and 
a Jordan matrix Zx> which has only zeros in the main diagonal, (accordingly to 
PQ, chapter 8.7) So that X is isomorphic to a quaternion then Zx is of type 
or lfT2|). In the first case X is the Zero matrix or isomorphic with it X is the 
Zero quaternion (the trivial solution). 

/ 1 \ 

In the second case we have X = Tj^X -1 = X 1 

\ j 

was computed, then we can compute a special matrix according to 1(13)1 . ((14)1 or 
(fT5l). Hence the solution is 



If the vector x 



X = { 



x = 0, 



= 0} 



(31) 



f x ° \ 

Xl 
V *3 J 

The trivial solution is included in this solution set. In contrast to the solution 
set at non-singular quaternions we do not have f spheres but only one. 
2. a ^ and a 2 ^ 0. (Case ®) 

So that 21 is singular, we have the condition a 2 , + a 2 = 0. Then 21 has the 
eigenvalues 2ao and and the Jordan decomposition is 21 = ilJU -1 with Z = 
/ 2a Q \ 
2a 




\ / 



If we restrict to solution matrizes that are isomorphic to a 



quaternion, this leads to type J5]) and we get: X = 
Here we have n different values for | -\/2oo. 



2 a a 



Example: Let A = 



X, = 



-8 + 8i \ 
-16 + 8i 
-8 + 24i 



V 


24 4 


-16i J 


1 1H 




( 


1 + 


2i 




3H 


- i 


x 2 = i 

- 2 


V 2- 


3i J 


V 



n = 3 The solutions are: 



/ -1- V3+(-l + V3> \ 
-l-2V3+(-2 + \/3)« 
-3- Vl+(-l + 3V3)i 



/ -l + V3+(-i-V3)t \ 
-l + 2V3 + (-2- V3)i 
-3 + V3~+ (-1 -3\/3)i 

V -l-3V3 + (3-2V3> / 



3. a ^ and a 2 = 0. (Kategorie (H))) 

So that 21 is singular then is a = 0. It follows the jordan normal form of 21. 
/ 1 \ 

"" and it according to CE3]), O or (ffU)- a is 



21 = itfiT 1 mit Z 



1 
\ / 
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nilpotent of index 2. The solution matrix X has to be nilpotent of index x with 
n < x < 2n. We have n > 2, therefore x > 3. It follows that X cannot be a 
quaternion isomorphic matrix. The equation X" = A is insoluble. Example: 
/ \ 




5 Summary 

Now the results are displaied shortly: The quaternion equation X™ = A, n > 2, 
integer, has the the following solution sets: 

(l\ 

, ao + 



la. A = a E 



a 






a) X = n different values of ^/clq are possible. 

b) Let wi and u>2 be 2 different n-th roots of ao (w™ = w% = a but w\ ^ w^)- 
In order to compute X we have to consider all combinations of w\ and w 2 ■ 

'' " ' ! 1 x = 5(tt>i + w 2 ), x 2 = x{ + x\ + a;§ = -|(tui - w 2 ) 2 }; (2) 



X = { 



Kombinationen 
lb. A = a E = a 



x = 0, 
2a. A = 



x = x 

«o 
a 



/ 1 \ 



V oj 

x 2 



a = 







, mit a ^ 0, a 2 7^ und a 2 + a 2 , 7^ 0. 



X = 



a 



2j ' 



a 



1 
2 

j__ 

'2i v 



with a 7^ 0, a 2 7^ and a 2 + a§ = (singular case) 



According to the multiplicity of the n-th root there are n 2 solutions to X. 

fl ° ^ . with a 0. r 2 -J- " — < - 2 ' - 2 
a 

2a / \ a 

According to the multiplicity of the n-th root there are n solutions to X. 
a a ° ] . with a - 0. f >-' - 




3a. A = 



, with a 7^ 0, a 2 = and ao 7^ 



J n V°o 

According to the multiplicity of the n-th root there are n solutions to X. 



3b. A 



ao 



, with a 7^ 0, a =0 and ao = (singular case) 



In this case we do not have a solution. 
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